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A low CMB variance in the WMAP data 
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ABSTRACT 

We have estimated the CMB variance from the three-year WMAP data, finding a 
value which is significantly lower than the one expected from Gaussian simulations us- 
ing the WMAP best-fit cosmological model, at a significance level of 98.7 per cent. This 
result is even more prominent if we consider only the north ecliptic hemisphere (99.8 
per cent). Different analyses have been performed in order to identify a possible origin 
for this anomaly. In particular we have studied the behaviour of single radiometers 
and single year data as well as the effect of residual foregrounds and 1/f noise, finding 
that none of these possibilities can explain the low value of the variance. We have also 
tested the effect of varying the cosmological parameters, finding that the estimated 
CMB variance tends to favour higher values of n s than the one of the WMAP best-fit 
model. In addition, we have also tested the consistency between the estimated CMB 
variance and the actual measured CMB power spectrum of the WMAP data, finding 
a strong discrepancy. A possible interpretation of this result could be a deviation from 
Gaussianity and/or isotropy of the CMB. 
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1 INTRODUCTION 

The study of the anisotropies of the cosmic microwave back- 
ground (CMB) constitutes one of the most powerful tools 
of cosmology, providing us with very valuable information 
about the origin and evolution of the universe. In particu- 
lar, the determination of the CMB power spectrum allows 
one to put tight constraints on the cosmological parameters. 
In addition, the study of the CMB temperature distribution 
provides us with a powerful test of the standard inflation- 
ary theory, since this predicts that the CMB fluctuations 
should follow a homogeneous and isotropic Gaussian distri- 
bution w hereas alternative th eories - such as non-st andard 
inflat ion l|Bartolo et al.ll2~004l ) or topological defects (|Durred 
1 19991 ) - give rise to non-Gaussian fluctuations. 

Given the interest of the subject, a large num- 
ber of techniques have been proposed for study- 
ing the Gaussianity and isotropy of the CMB in- 
cludi ng, among others, th e Mi nkowski function- 



timators ([Ba rreiro, Martinez-Gonzalez fc Sanzl 



[ Pore. Colombi fc Bouchet] 
20061) extrema correlation function 
1 19951 ; iBarreiro et al.l 1998; Heavens 



20031; iMonteserin et all 



2001 : 
20051 



dNaselskv fc Novikovl 
Shethl | l999h. wavelet 
analysis (iHobson. Jones fc Lasenbvl 1 19991: IBarreiro et al 



I 200C : IBarreiro fc Hobsonll200ll : l Aghanim. Forni fc Bouchet 
12001 ) , bipolar power spectrum llHaiian fc Souradeedl2003[ r 
phase analysis (Chiang. Nas elskv fc Co les 2004) and good- 
ness of fit tests ( Aliaga et al.l 120051 ; iRubino-Martm et al.l 
l200fj ; ICurto et aljkoOTt ) 



als dColesl Il988l; iGott et al.1 1 1990h . the bispectrum 
1 Ferreira. Magueiio fc Go rski 1998; Heavens 1998; Magueiiol 



2000h . properties of hot and co l d spo ts ( Coles fc Barrow! 
19871 : iMarti'nez-Gonzalez et all |2000| ) . geometrical 



Many non- Gaussianity analyses have been performed 
using the best CMB data available up to date, which 
have been provided by the NASA Wilkinson Microwave 
Anisotropy Probe (WMAP) satellite l|Bennett et alj|2003al . 
ISpergel et al] 120071 ). The WMAP team carried out a study 
of the CMB temperat ure distribution of the first year data 
l|Komatsu et alj|2003l ) finding that the data were consistent 
with Gaussianity. However, in subsequent works, a number 
of unexpected results regarding the Gaussianity and/or 
isotropy of the CMB were reported, including anomalies 
related to low multipoles dde Oliveira-Costa et all |2004|; 
ICopi et all 120041 : lLand fc Magueiiol l2005al ; ISchwarz et all 
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20041'). north/south asymmetrie s dEriksen et al.l 2004 



Hansen. Bandav fc G orski l2004al; lHansen et al.1 1 2004b : 



Eriksen et al.l 120051 : lLand fc Magueiid l2005bl : iBernui et ail 



a cold spot in the southern hemisphere rtVielva et al 



alignment Jwiaux et al 



Mukheriee fc Wand [iooi; ICruz et al 



McEwen et al.l 120051 : ICavon et al.1 I2005T ). structure 

20061). phase correlations 
120031 ; IColes et all 
hot and cold spots 

I Larson fc Wandeltll20o4 ). 

After the release of the three-year WMAP data 
March 2006, similar analyses have been carried 
confirming the presence of the anomalies in the 



dChiang. Naselskv fc Verkhodanovl 
I2004T) and to the am plitude of 



data ((Martinez-Gonzalez et al. 


20061; McEwen et al.l 


200£ 


Vielva et al.1 20071: Cruz et al. 


120071: CoDi et all 


2007 


Bernui et al.ll2007; Eriksen et al.ll2007; Wiaux et al.li2008). 



In this work, we report a new anomaly in the WMAP 
data: a significantly low value of the CMB variance with 
respect to the one expected for the WMAP best-fit model. 
The outline of the paper is as follows. In fj2]several quantities 
related to the one-point density function (1-pdf) are stud- 
ied for the WMAP data, including the dispersion, skewness, 
kurtosis and Kolmogorov-Smirnov distance. In Sj3]the CMB 
dispersion from the WMAP data is estimated, finding that 
its value is significantly low. The effect on our results of pos- 
sible residual foregrounds or systematics are studied in Sj4] In 
fj5] possible explanations for the anomaly are discussed, in- 
cluding the modification of the cosmological parameters and 
a deviation of the CMB from Gaussianity and/or isotropy. 
Finally, our conclusions are presented in §6] 



2 TESTS OF GAUSSIANITY BASED ON THE 
1-PDF 

The NASA WMAP satellite was launched in the summer of 
2001. The first-year and three-years results were presented 
in February 2003 and March 2006 respectively. WMAP ob- 
serves at five frequency bands: K (22.8 GHz, one receiver), 
Ka (33.0 GHz, one receiver), Q (40.7 GHz, two receivers), 
V (60.8 GHz, two receivers) and W (93.5 GHz, four re- 
ceivers). All the data and products generated by WMAP 
can be found at the Legacy Archive for Microwave Back- 
ground Data Analysis (LAMBDA ) web siteQ] 

Following the WMAP team i|Komatsu et al.ll2003l ) we 
have performed our non- Gaussianity analyses on a noise 
weighted average of the Q, V and W receivers, which pro- 
vides a CMB map where the signal-to-noise ratio has been 
increased. The K and Ka receivers are not included in the 
combined map since they are largely contaminated by Galac- 
tic foregrounds. In addition, the WMAP team has reduced 
the Galactic contamination present in the Q, V and W map 
by performing a for eground template fit as described in 
IHinshaw et~al] (|2007h . In particular, the comb ined Q+V+W 
map is constructed using these clean data as (jBennett et al.l 
l2003bl ): 



10 

T c (x) = ^T j (x) 

i=3 



w 



(1) 



where x gives the position in the sky and the index j cor- 
respond to the different receivers of the Q, V and W bands 
(i.e., the indices 3 to 10 refer to the Ql, Q2, VI, V2, Wl, 
W2, W3 and W4 radiometers respectively). The noise weight 



uij(x) is defined as: 



Wi x = 



£,=3%( X ) 



(2) 



where aoj is the n oise dispersion per o bservation for each ra- 
diometer given bv lJarosik et al.ll2007l and iVj(x) is the num- 
ber of observations made by the receiver j at the position in 
the sky x. 

Although the data are provided at a HEALPi>0 
|G6rski et al.1 [2005) resolution of n S ide = 512, we degrade 
the data down to n 3i d e = 256 since the smallest scales are 
dominated by noise. In addition, in order to avoid the strong 
contamination present at the Galactic plane and the emis- 
sion coming from extragalactic point sources, only t he data 
outside the WMAP KpO mask (IHinshaw et al.ll2007l ). which 
corresponds approximately to 76 per cent of the sky, has 
been used. Finally, the monopole and dipole outside the KpO 
mask have been removed. 

In order to apply our Gaussianity test, we construct 
the normalised temperature w(x), which is obtained by di- 
viding the data at each pixel by its corresponding expected 
dispersion. Since the contribution from residual foregrounds 
in the clean WMAP combined map outside the KpO is ex- 
pected to be very small, the main contribution to the data 
dispersion comes from the CMB signal and the instrumental 
noise. The WMAP noise is very well approximated by Gaus- 
sian white noise at each pixel characterised by a dispersion 
a n (x). Therefore, the normalised temperature is given by 



it(x) = 



T c fx) 



(3) 



^Jal +<jI (x) 

where <jq is the CMB dispersion^ and has been estimated as 



cro 



\ 



21 + 1 

i=2 



Cf 



(4) 



http: / /cmbdata. gsfc.nasa.gov 



The Ci's in the previous equation correspond to the power 
spectrum of the combined map assuming the best-fit modeQ 
to the WMAP data and taking into account the differ- 
ent beams and noise weights used in the combination (see 



2 http://healpix.jpl.nasa.gov 

3 Since the CMB is an homogeneous and isotropic field, cro is 
constant over the sky. However, strictly speaking, there is a small 
dependence of ctq with the position in the WMAP combined map 
due to the different beams and noise weights used in the combi- 
nation. Nonetheless, these differences are very small and one can 
assume, to a very good approximation, that the signal dispersion 
is constant over the sky. 

4 The CMB simulations have been performed using the best- 
fit model — for the case of a ACDM model using only WMAP 
data — given in the first v ersion of the paper s of the three- 
year WMAP data realease llSpereel et al.1 120071, see table 5 of 
arXiv:astro-ph/0603449vl). Although the values of the cosmo- 
logical parameters have been slightly modified in the final version 
of the WMAP papers, we have tested that these small changes 
do not affect the results of this work. 
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Table 1. The dispersion, skewness, kurtosis and Kolmogorov- 
Smirnov distance dKS are given for the normalised WMAP data 
as well as for the northern and southern ecliptic hemispheres (in 
all cases considering only pixels outside the KpO mask). The sig- 
nificance of each quantity has been calculated from 1000 Gaussian 
realisations as the percentage of simulations with a value larger 
or equal to that obtained from the data. 



Estimator 


WMAP value 


Signif. (%) 


All pixels outside KpO mask 


dispersion 


9.55 x lO" 1 


97.8 


skewness 


-1.91 x 10~ 2 


74.7 


kurtosis 


1.49 x 10" 2 


20.9 


dKS 


1.21 x 10~ 2 


7.0 


Northern ecliptic hemisphere 


dispersion 


9.35 x 10" 1 


99.4 


skewness 


2.17 x 10" 2 


28.1 


kurtosis 


2.18 x 10" 2 


20.1 


dKS 


1.89 x 10~ 2 


25.2 


Southern ecliptic hemisphere 


dispersion 


9.76 x lO" 1 


71.8 


skewness 


-5.85 x 10~ 2 


92.8 


kurtosis 


7.08 x 10- 4 


35.4 


dKS 


1.05 x 10" 2 


67.3 



Appendix [XJ. The maximum t that we have considered is 

If the data follow a Gaussian distribution - and our as- 
sumptions about the underlying cosmological model and the 
instrumental noise are correct - the normalised temperature 
u(x) should follow a Gaussian distribution with zero mean 
and unit dispersion N(0,1). 

We have calculated different estimators for the nor- 
malised data outside the KpO mask: dispersion, skewness 
and kurtosis. In addition, the data have been compared to a 
N(0,1) distribution through a Kolmogorov-Smirnov test. In 
order to assign a significance to these quantities, the same 
analysis has been applied to 1000 Gaussian simulations of 
the WMAP combined data. To reproduce the data as closely 
as possible, the simulations have been generated at resolu- 
tion n s ide=512 for each radiometer and then processed in 
the same way as the data to obtain the simulated combined 
maps. Results are shown in Table [T] We find that the skew- 
ness and kurtosis are consistent with the results obtained 
from Gaussian simulations. However, the dispersion of the 
normalised data is significantly lower than expected, at the 
level of 97.8 per cent (i.e., only 22 out of the 1000 simula- 
tions have lower values than the one found for the data). 
Regarding the Kolmogorov-Smirnov distance (dies), we find 
that only 7 per cent of the simulations have a larger value 
of dKS than the data. This could be related to the fact that 
the data dispersion is lower than expected, what would tend 
to produce a deviation of the normalised temperature from 
a N(0,1) distribution. These deviations can also be seen in 
the top panel of Figure [TJ which compares the histogram 
obtained from the normalised WMAP data with the one 
obtained averaging over 1000 Gaussian simulations. A de- 
viation of the data histogram (dashed line) from the one 



obtained with simulations (solid line) is apparent due to the 
small value of the dispersion of the normalised data. We 
would like to mention that the same analysis has also been 
performed on the first year WMAP data, finding again a 
low value of the dispersion (at a significance level of 99.5 
per cent). 

We may wonder if errors in the estimation of the 
WMAP noise level could affect the results of table [1] In 
particular, taking into account equation ©, an overestima- 
tion of the noise level would tend to bias the dispersion of 
the normalised temp erature towards values lower than unity. 
Jarosik et all (|2007h use two different methods to estimate 
the noise level of each WMAP radiometer, finding that both 
methods agree within 0.3 per cent. This value can be con- 
sidered as an indication of the level of error expected in 
the estimation of the WMAP noise. Even if we consider an 
overestimation of 1 per cent in the noise dispersion and re- 
calculate, accordingly, the normalised temperature for the 
data, we still find a significant deviation for the dispersion. 
In particular, the new value of the dispersion is 0.956, which 
corresponds to a significance of 97.4 per cent. Therefore, the 
level of errors expected in the estimation of the noise level 
cannot explain this anomaly. 

In order to investigate further the origin of this devia- 
tion, we have repeated the same analysis in the northern and 
southern ecliptic hemispheres, considering only those pixels 
outside the KpO mask (which correspond to a sky fraction of 
38.8 and 37.0 per cent for the northern and southern hemi- 
spheres, respectively). Results are shown in Table [T] The 
values of skewness, kurtosis and dKS remain compatible with 
Gaussian simulations. However, we find an asymmetry in the 
behaviour of the dispersion: while the data dispersion in the 
southern ecliptic hemisphere is compatible with the simula- 
tions (at the 71.8 per cent level), a clear deviation is found 
in the north, at the 99.4 per cent significance level (i.e., only 
6 out of 1000 simulations had lower dispersion values). The 
middle (bottom) panel of Figure [T] shows the histogram of 
the data in the northern (southern) ecliptic hemisphere com- 
pared to the averaged histogram obtained from Gaussian 
simulations in the same sky region. The histogram obtained 
from the northern data shows again a systematic deviation 
from the average one, which is even larger than in the case of 
the whole sky, due to the lower value of the dispersion found 
for the data of the northern ecliptic hemisphere. Conversely, 
the histogram obtained from the southern data is closer to 
the average value of the Gaussian simulations. 



3 STUDY OF THE CMB DISPERSION ao 

Given the deviation found in the dispersion of the nor- 
malised temperature, we have performed a more detailed 
study of the CMB dispersion ao. Let us recall that the dis- 
persion of the data cr(x) has contributions from both the 
CMB signal and the instrumental noise: 

cr(x) = v/a 2 + a 2 (x) (5) 

We have estimated the value of ao by allowing this quantity 
to vary in equation <(3j and compare the different obtained 
normalised temperatures to a Gaussian of zero mean and 
unit dispersion with the Kolmogorov-Smirnov test. Our es- 
timation of the CMB dispersion is given by the value of ao 
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Figure 1. The top panel shows the histogram of the normalised 
WMAP data outside the KpO mask (dashed red line) compared 
to the averaged histogram obtained from 1000 simulations (solid 
black line). The error bars indicate the dispersion obtained from 
simulations. The middle and bottom panels give the same his- 
tograms for the northern and southern ecliptic hemispheres (only 
considering pixels outside the KpO) respectively. 



0.1 




0.08 0.085 0.09 0.095 0.1 

true a (mK) 

Figure 2. Estimated versus true dispersion (in mK) obtained 
from 1000 simulations of the WMAP combined map outside the 
KpO mask. For illustration, the best fit to a straight line is also 
shown (solid red line), which corresponds to ctq = 1.027cro — 0.002. 



that produces the minimum cIks ■ We have tested with simu- 
lations that this estimator performs well, as it will be shown 
below for the different considered cases. 

Using this method, we have estimated a value of ao = 
8.16 x 10" 2 mK for the WMAP combined map. The bias and 
error of this estimation has been tested using 1000 CMB sim- 
ulations. Figure [2] shows the true CMB dispersion of each 
simulation outside the KpO mask versus the one estimated 
with our method from the simulated data in the same re- 
gion of the sky. For illustration, the best fit to a straight line 
is also plotted, showing a very good correlation between the 
true and estimated dispersions. The bias (obtained as the av- 
erage difference of the true minus the estimated dispersion) 
and error (obtained as the dispersion of the same difference) 
of the estimator are b = —6.9 x 10 -5 mK and e = 4.2 x 10~ 4 
mK respectively, i.e. the error of our estimator is at the level 
of ~0.5 per cent. In addition, the bias is also significantly 
smaller than the typical value expected for the CMB disper- 
sion and therefore our estimator can be considered unbiased 
to a very good approximation. An additional bias in the es- 
timated dispersion could also be present if the assumed level 
of instrumental noise were incorrect. In particular, an over- 
estimation of the noise would imply an underestimation of 
the CMB dispersion. To check this possibility, we have also 
estimated cr as (TiT 2 ), with 7\ = Ql + VI + Wl + W2 
and T 2 = Ql + V2 + W3 + W4, finding a value which is 
consistent with the previous estimation within the la error. 
Note that since instrumental noise is uncorrelated between 
channels, this latter estimation is not affected by a possible 
bias in the noise level. 

For comparison, we have also calculated the mean ero 
from 20000 noiseless CMB Gaussian simulations of the 
WMAP combined data, finding a value of 8.60 x 10~ 2 mK, 
clearly higher than the dispersion estimated for the data. 
In order to quantify the significance of this result, we have 
also obtained the distribution of variances from the same 
simulations. We find that 98.7 per cent of the simulations 
have a larger value of the variance (or, equivalently, of the 
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Table 2. The second column gives the CMB dispersion (in mK) estimated for the WMAP combined data obtained using pixels of three 
different regions of the sky: outside the KpO, the northern ecliptic hemisphere outside the KpO and the southern ecliptic hemisphere 
outside the KpO. In the third column, the mean dispersion value of the model, obtained averaging over 20000 Gaussian simulations, is 
given. Significances, that have been calculated as the percentage of simulations with values larger than the one found for the data using 
20000 Gaussian CMB simulations, are given in the fourth column. Finally, in the last two columns, the bias (obtained as the average 
differences of the true minus the estimated dispersion) and error of the estimator, both in mK, are given, which have been obtained from 
1000 Gaussian simulations. 



Region 






model 
"0 




Signif. 


Bias 




Error 




KpO 


8.16 x 10 


-2 


8.60 x 10" 


2 


98.7 


-6.9 x 10- 


5 


4.2 X 10" 


-4 


north 


7.97 x 10 


-2 


8.59 x 10" 


-2 


99.8 


-1.6 x 10" 


-4 


6.8 x 10" 


-4 


south 


8.37 x 10 


-2 


8.59 x 10" 


-2 


75.7 


-2.0 x 10" 


-4 


7.2 x 10" 


-4 



dispersion) than the one found for the WMAP data. These 
results are summarised in Table [2] 

It is interesting to point out that, for a Gaussian CMB, 
the distribution of the correlation function C( 8) for a given 
powe r spectrum can be analytically calculated (|Cav6n et al.l 
Il99ll ). In particular, we are interested in obtaining the theo- 
retical distribution of t he variance, i.e., the correlation func- 
tion at 8 = 0. Following ICavon et alj (|l99ll ). the cumulative 
function of the variance can be calculated as 



n[ 

1=0 

sin [B (i) + ta% 



B ® = E 



t 

21 + 1 



dt 



arctan (— 2t<r^) 



(6) 



where ae — ^ contains the power spectrum dependence. 
Calculating the derivative of the cumulative function with 
respect to Oq, we can also obtain the probability distribution 
of the variance: 



( *\ 1 
p{ff ) = - 

7T 



tmax 

n[< 



1 + 4t 2 ai) 



[B (t) + tal] dt(7) 



The previous distributions are valid when considering the 
variance of the whole sky. However, the presence of the KpO 
mask restricts the fraction of the sky (~ 76 per cent) used in 
the analysis. Since obtaining the significance of a given esti- 
mated CMB variance from simulations is significantly slower 
than calculating it from the cumulative function (0 , we have 
tested whether the theoretical distributions are a good ap- 
proximation in our case. The left panel of Figure [3] shows 
the cumulative function of erg (black solid line) obtained us- 
ing the theoretical power spectrum of the WMAP combined 
map calculated in AppendixfS] The dotted-dashed blue line 
corresponds to the distribution of the variance obtained from 
60000 Gaussian simulations over the whole sky. As expected, 
the agreement in this case is very good. Finally, the red 
dashed line gives the distribution of the variance obtained 
from the same simulations considering only those pixels out- 
side the KpO mask. We can appreciate some small differences 
between this curve and the one obtained from equation (|6]), 
due to the fact that we are applying a mask, but the general 
agreement is very good. The maximum difference between 
both cumulative functions is found in the central part of 
the plot and is less than 0.04, whereas the differences in 



the tails of the distributions are even smaller. Therefore, 
the theoretical cumulative function provides a simple and 
good approximation to obtain the significance of the esti- 
mated CMB variance in the considered case. In particular, 
using the theoretical distribution we obtain a significance of 
99.3 per cent for the estimated 00, that should be compared 
with the value of 98.7 obtained from simulations. Similar 
conclusions can be derived from the right panel of Figure |3j 
that shows the probability distribution of the variance for 
the same cases as before: theoretical (solid black line), from 
simulations using the whole sky (blue dotted-dashed line) 
and from simulations considering pixels outside KpO (red 
dashed-line) . 

Taking into account the results found in the previ- 
ous section, we have also estimated the CMB dispersion in 
the northern and southern ecliptic hemispheres (considering 
only those pixels outside the KpO mask), obtaining values 
of 7.97xl0~ 2 and 8.37xl0 -2 mK, respectively. The corre- 
sponding biases, errors, mean dispersion of the model and 
significances are given in Table [2] As expected, the bias and 
error of the estimator has increased slightly with respect to 
the ones obtained using all pixels outside the KpO mask, due 
to the smaller fraction of sky considered, but they are still 
small (less than 1 per cent with respect to <xo). In addition, 
it should be noted that this bias is negative, implying that, 
if anything, the CMB dispersion would tend to be overesti- 
mated and, therefore, the true 00 would be even lower. 

We remark that the variance of the northern hemisphere 
is significantly low (only 0.2 per cent of the Gaussian simula- 
tions have lower values) whereas the southern hemisphere is 
compatible with the simulations. This is consistent with the 
results found in §[21 which seem to indicate that the anomaly 
found in the WMAP variance comes from the northern eclip- 
tic hemisphere. This result is also in agreement with the 
lack of power found in the northern e cliptic hemisphere of 
the WMAP data by previous works (lEriksen et al. 1 12004 
lHansen et al.ll2004al : Eriksen et alj|2005l ). 

We have also studied the ratio of <r 2 in the northern 
ecliptic hemisphere over <tq in the southern one for the 
WMAP data, finding that this quantity is consistent with 
what is expected from Gaussian simulations (with a signif- 
icance of 88.7 per cent obtained from 20000 simulations). 
This can be explained because, although the value of 00 in 
the southern hemisphere is consistent with Gaussianity, this 
quantity is lower than the average dispersion expected from 
simulations, and this fact makes the north-south ratio com- 
patible with the Gaussian model. 
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Figure 3. Left panel: theoretical cumulative function of the CMB variance (solid black line), compared to the same functions obtained 
from 60000 Gaussian simulations over the whole sky (blue dotted-dashed line) and using only the pixels outside the KpO mask (red 
dashed line). Right panel: the probability distribution of the variance is given for the three same cases. For illustration, the solid vertical 
line in both panels indicates the value obtained from the WMAP data. The units of the variance are (mK) 2 . 



Finally, we have studied the effect that the low CMB 
quadrup ole measured by COBE dBennett et al.lll996f) and 
WMAP l|Bennett et al.ll2003al : Iffinshaw et al.l 120071 ) has in 
our results. Taking into account equation Q and the fact 
that the Ci 's are roughly proportional to l/£(£+l), it is easy 
to see that the largest contributions to the CMB variance 
come from the lower multipoles and, in particular, from the 
quadrupole. In order to test if the low CMB quadrupole can 
be the origin of the low CMB variance we have repeated 
our analysis after removing the best-fit quadrupole (outside 
the KpO mask) from the data. With this procedure we find 
(To = 8.07 x 10 -2 mK which corresponds to a theoretical 
significance of 96.9 per cent (obtained using equation ((6} 
setting C2 =0). In addition, we have obtained the signif- 
icance from 20000 simulations, where the quadrupole has 
been removed, finding a significance of 95.0 per cent. The 
same analysis has also been performed in the northern and 
southern ecliptic hemispheres, estimating for the dispersion 
7.91 x 10" 2 mK and 8.30 x 10" 2 mK, respectively. The signif- 
icances of these values, obtained from 20000 Gaussian sim- 
ulations, are 99.2 and 44.3 per cent, respectively. Therefore, 
after removing the quadrupole, although the significances 
are slightly reduced, the CMB variance of the whole sky 
and that of the northern ecliptic hemisphere are still anoma- 
lously low. Thus, the low quadrupole cannot explain by itself 
these deviations. 



4 DISCUSSION 

We may wonder if the anomalous value found for the CMB 
dispersion in the WMAP data may have an extrinsic ori- 
gin, i.e., to be due to the presence of some systematics or 
foreground residuals. In order to clarify this point, we have 
performed some further tests including an analysis of each 
individual radiometer ( £|4. ip as well as the study of the effect 
of 1/f noise ( i|4.2|l and foreground residuals ( i|4.3l) . 



4.1 Single radiometer analysis 

In order to check if the low value of ero found in the WMAP 
data could be explained by the presence of an anomalous ra- 
diometer, we have estimated the CMB dispersion from the 
maps of each single radiometer following the method ex- 
plained in fJ3] The results are summarised in Table The 
mean dispersion a ™ odel and the significances (see Sj3} have 
been obtained analytic ally using the power spectrum of the 
WMAP best-fit model (|Spergel et al.ll2007l ) and taking into 
account the different beam window functions of each ra- 
diometer as well as the window pixel function of the consid- 
ered resolution (N si< ie ~ 256). The biases and errors of the 
estimator have been obtained from 100 Gaussian simulations 
and are at a level less than 1 per cent in all cases. 

As seen in Table [3] a different value of 00 has been ob- 
tained for each radiometer, which is expected due to the 
different beam window functions of each radiometer as well 
as to the statistical error of our estimator. In order to check 
if the dispersion of each radiometer is consistent with each 
other, and that no single radiometer has an anomalous dis- 
persion that could be responsible for the low value of <to 
found in the combined map, we have performed the follow- 
ing test. First, we have obtained the mean, d, and dispersion, 
error(d), of the difference between the dispersion of the com- 
bined map minus that of each single radiometer from 100 
simulations, where the map dispersions were obtained using 
our estimator. We have then added the corresponding mean 
difference d to the &o of each radiometer, which provides an 
estimation of <to for the combined map (given in the sixth 
column of Table |3J). If all the radiometer maps are consistent 
with each other, all these quantities should be compatible 
among them within the range allowed by the dispersion of 
the calculated differences (given in the seventh column of 
Table [3} as well as with the 60 derived for the WMAP com- 
bined data. From Table [3] we see that consistency is found 
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Table 3. CMB dispersion (in mK) estimated from maps of individual radiometers outside the KpO mask. The mean dispersion of the 
model and the significances have been obtained analytically. The corresponding bias and error (in mK) have been obtained from 100 
Gaussian simulations. For comparison, the same results obtained for the WMAP combined map are also shown. An estimation of the 
(To of the combined map from each of the radiometers and its corresponding error are also given in the last two columns (see text for 
details). 



Radiometer 


O"0 




model 




Signif. (%) 


Bias 




Error 




o-o + d 




Error (d) 


Ql 


7.93 x 10- 


-2 


8.34 x 10- 


-2 


98.0 


-4.9 X 10 


-5 


4.1 X 10" 


-4 


8.21 x 10" 


-2 


3.7 x 10" 4 


Q2 


7.88 x 10" 


-2 


8.29 x 10- 


-2 


97.9 


4.1 x 10- 


6 


4.5 x 10" 


-4 


8.21 x 10- 


-2 


3.7 x 10- 4 


VI 


8.42 x 10- 


-2 


8.84 x 10" 


-2 


98.9 


1.0 x 10- 


5 


4.4 x 10- 


-4 


8.19 x 10- 


-2 


3.7 x 10~ 4 


V2 


8.46 x 10- 


-2 


8.88 x 10" 


-2 


98.8 


-8.4 x 10 


-6 


3.7 x 10- 


-4 


8.19 x 10- 


-2 


3.6 x 10" 4 


Wl 


8.66 x 10" 


-2 


9.05 x 10" 


-2 


98.3 


-2.5 x 10 


-5 


5.4 x 10" 


-4 


8.20 x 10- 


-2 


5.7 x 10" 4 


W2 


8.48 x 10- 


-2 


8.87 x 10" 


-2 


98.1 


-8.7 x 10 


-6 


5.8 x 10" 


-4 


8.21 x 10" 


-2 


6.6 x 10- 4 


W3 


8.47 x 10- 


-2 


8.90 x 10- 


-2 


99.1 


1.1 x 10- 


4 


6.8 x 10" 


-4 


8.16 x 10- 


-2 


6.5 x 10~ 4 


W4 


8.74 x 10- 


-2 


9.07 x 10" 


-2 


95.9 


8.2 x 10 - 


5 


5.8 x 10- 


-4 


8.26 x 10- 


-2 


5.9 x 1Q- 4 


Combined 


8.16 x 10- 


-2 


8.62 x 10- 


-2 


99.3 


-6.9 x 10 


-5 


4.2 x 10- 


-4 









and that no single radiometer has an anomalous dispersion 
Moreover, the significances obtained for the estimated dis- 
persion for all the radiometers are significantly high, con- 
firming that the anomaly cannot be attributed to a single 
radiometer. 

In addition, we have repeated the same analysis for 
all the radiometers for each individual year of data. We 
find that the estimated CMB dispersion is reasonably stable 
along the three years of data for all the radiometers as well 
as for the combined map. This shows that the low ao cannot 
be explained as the result of some anomaly in a single year 
of data. 



4.2 1/f noise 

In our analysis, it has been assumed that the WMAP in- 
strumental noise is a Gaussian white noise characterised by 
a dispersion er„(x) which varies with the position of the 
sky. Although this seems to be a very good approxima- 
tion, the WMAP data also contain a small level of 1/f noise 
|jarosik et alj|2003t [20071 ) and therefore we may wonder if 
this could be affecting our estimation of ao- 

In order to test the effect of 1/f noise we need realistic 
noise simulations that reproduce the WMAP data pipeline. 
As far as we know, this type of simulations are not available 

5 Although consistency is found, it is observed that the disper- 
sions derived from the radiometers are in general higher than 
that derived directly from the WMAP combined map. This small 
bias could be due to the presence of residual foregrounds and 
1/f noise in the single radiometer maps, whose effect, as will be 
shown in sections i|4,2l and £|4.3I is to produce an overestimation 
of the dispersion in these maps. This is consistent with the fact 
that the radiometers that have a dispersion that deviates more 
from the one of the combined map (although always within the 
la error) are Ql and Q2, which are expected to be the channels 
more affected by residuals foregrounds, and W4, which present s 
the largest 1/f noise of all the radiometers l ljarosik et al.|[2003l) . 
These possible residuals are expected to be reduced, relatively to 
the CMB signal, in the combined map and, therefore, the direct 
estimation of the ao from the combined data would tend to give 
a lower (less biased) value than that obtained from each of the 
radiometers. 



for the three-year data but, however, the WMAP team did 
provide 110 of these simulations for each radiometer for the 
first year of data that include 1 /f noise as well as all known 
radiometric effects (for details see the information in the 
LAMBDA web site). 

Therefore, we have carried out the following test. First 
of all, we have combined the simulations of the different ra- 
diometers in the appropriate way in order to obtain 110 re- 
alistic noise simulations of the first-year WMAP combined 
map. We have then obtained the realistic noise dispersion 
map from these simulations. Since the number of simula- 
tions is relatively small, if we try to obtain the dispersion at 
each pixel, this map will be very noisy. Therefore, we have 
assumed that this quantity varies smoothly over the map 
and that neighbouring pixels would have the same level of 
noise. In particular, we have obtained the dispersion from 
16 neighbouring pixels along the 110 simulations. The mean 
value of the ratio of this dispersion map over the one con- 
structed assuming that only white noise is present is 1.02, 
which is an indication of the small level of 1/f noise present 
in the data. 

The next step is to construct a new set of 110 realistic 
simulations by adding the combined realistic noise simula- 
tions to Gaussian CMB realizations of the WMAP combined 
first-year data. We have then applied our dispersion estima- 
tor to this new set of simulations and study its bias and error 
in 2 different cases: 1) including in equation (J5]) the value of 
cr n (x) obtained from the realistic noise simulations and 2) 
assuming that we have only white noise. Interestingly, the 
errors in both cases are very similar (~4x 10 -4 mK), but 
the biases are significantly different. In the first case, when 
the noise is correctly characterised, the bias is very small 
(—3.4 x 10~ 5 mK), however, when only white noise is erro- 
neously assumed, the bias increases to a value of —4.6 x 10~ 4 
mK. This means that, if 1/f noise were present in the data 
and we did not properly take it into account, our estimation 
of the CMB dispersion would tend to shift towards higher 
values and, therefore, this effect could not explain an anoma- 
lous low value of ao- This can be understood, since, by not 
including the 1/f noise, we are underestimating the level of 
instrumental noise and this extra power would tend to be 
compensated by increasing the signal dispersion. As a fur- 
ther test, we have estimated the CMB dispersion from the 
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WMAP first-year data using the realistic noise dispersion 
map, finding a value of 8.12 x 10 -2 mK to be compared 
with the larger dispersion obtained assuming white noise 
(8.17 x 1(T 2 mK). 

Since there are not realistic noise simulations available 
for the WMAP three-year data, we cannot repeat this test 
for the complete set of data. However, we can qualitatively 
check that the same argument is valid for the three-year 
data using the reduced resolution inverse noise covariance 
matrices provided in the LAMBDA web site for the one and 
three-year data. The diagonal elements of the (direct) noise 
covariance matrices are dominated by white noise, whereas 
the off-diagonal elements come from 1/f noise. We have ob- 
tained the ratio of the diagonal and off-diagonal terms of 
the first-year data over the corresponding elements of the 
three- year data finding that, in both cases, this ratio peaks 
around 3. This means that the relative importance of the 1/f 
noise with respect to the white noise is approximately the 
same in both sets of data and, therefore, the effect that the 
1/f noise has in the three-year data should be very similar 
to the one found for the first-year data. 

Thus, residual 1/f noise cannot explain the low CMB 
variance found for the three-year WMAP data, since its un- 
accounted presence, if anything, would lead to higher values 

Of (TQ. 



4.3 Foreground residuals 

Although we are using the foreground cleaned maps pro- 
vided by the WMAP team and applying a Galactic plus 
point source mask, some foreground residuals may still be 
present in the data. Therefore, it is interesting to test if this 
possible contamination can affect our estimation of the CMB 
dispersion. 

An indication of a significant effect from foregrounds 
would be that the anomaly found in ao was frequency de- 
pendent. However, from Table [3] we can see that this is not 
the case, since ao is consistently low in the whole range of 
frequencies probed by Q, V and W (from 41 to 94 GHz). 

In order to study further the effect of foreground con- 
tamination, we have estimated the value of <to in the com- 
bined and individual radiometer maps without the fore- 
ground correction. The most relevant change is found for 
the Ql and Q2 maps (lowering the significances of ao down 
to 90.7 and 86.2 per cent respectively), whereas for the ra- 
diometers of the V and W bands, the changes are smaller. 
This is consistent with the fact that foreground emission 
is higher in the Q band and therefore the cleaning of the 
Ql and Q2 maps is important to improve the estimation of 
the CMB dispersion. The presence of residuals also trans- 
lates into a larger value for the ao of the combined map 
(8.23 x 1(T 2 mK). 

As one would expect, these results indicate that if fore- 
ground residuals were present in the data, this would tend 
to bias ao towards higher values. Therefore, the presence of 
foreground contamination cannot explain the anomaly found 
for the CMB dispersion. 



Table 4. Ratio between the maximum value of the likelihood, 
when varying each parameter within a range given by 2 times the 
errors of the WMAP best fit model, over the likelihood of the 
best fit model (column 2) and significance of the variance for the 
model correspoding to this maximum likelihood value (column 
3). 



Parameters 


Ratio 


Signif.(%) 


n b h 2 


2.1 


98.0 


Cl c h 2 


6.8 


90.3 


h 


3.6 


96.5 


T 


10.1 


68.5 


n s 


11.9 


79.0 


A 


14.1 


57.5 



5 POSSIBLE EXPLANATIONS OF THE 
ANOMALY 

In all the analyses carried out in this work, we have as- 
sumed that the CMB fluctuations are isotropic and follow a 
Gaussian distribution characterised by a given power spec- 
trum (which in turn is described by a particular cosmological 
model) . However, since we have found an anomaly regarding 
the variance that cannot be easily explained by foregrounds 
or systematics, we may wonder if some of these assump- 
tions are incorrect. In H5.1l we discuss if a different choice of 
cosmological parameters can make the estimated variance 
consistent with the model. In i]5.2l we study if this anomaly 
can be due to a more fundamental reason: a deviation of the 
CMB from Gaussianity and/or isotropy. 

5.1 Effect of cosmological parameters 

During this work, we have used for the CMB power spec- 
trum the best- fit model to the WM AP data assuming a flat 
ACDM model (|Spergel et al.ll2007h . However, we would like 
to test if by varying the values of the cosmological parame- 
ters - within the current uncertainties - the estimated CMB 
variance could be more consistent with the theoretical ex- 
pectations. 

In order to test this possibility, we have considered 
the effect of varying the following cosmological parame- 
ters: the physical baryon density Qbh 2 , the physical cold 
dark matter density Q. c h 2 , the reduced Hubble constant 
h, the reionization optical depth r, the spectral index n s 
and the amplitude of the density fluctuations A (with the 
constraint of having a flat ACDM model). Given that we 
have only one observable, the CMB variance oq, modify- 
ing all the cosmological parameters simultaneously would 
lead to strong degeneracies between them. Therefore, as a 
first approach, we have studied the effect of varying one 
single parameter (within a range given by 2 times the er- 
rors of the WMAP best-fit model) while leaving the others 
unchanged. In particular, we have constructed the corre- 
sponding one-dimensional likelihoods in the following way: 
for each set of par ameters, the CMB spectrum was calcu- 
lated using CAMB ([Lewis. Challinor fc LasenbylfeoOOT ) ; this 
power spectrum was then inserted in equation l)A8j) in or- 
der to get the Ce's of the WMAP combined map; finally, 
the probability of having the observed value of a 2 for that 
power spectrum was obtained using equation J7} . 
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We find that the estimated CMB variance favours lower 
values of flbh 2 , h and A than the ones of the best-fit model 
whereas higher values of Q c h 2 , n s and r are preferred. How- 
ever, not all parameters are equally sensitive to the variance 
as can be seen in Table [4] This table gives the ratio between 
the maximum value of the likelihood, when varying one pa- 
rameter in the considered range, over the likelihood of the 
WMAP best-fit model. The higher this ratio, the more sen- 
sitive the parameter to the CMB variance. We also show in 
Table [4] the significance of the variance for the model corre- 
spoding to this maximum likelihood value. We see that the 
compatibility of oa with the model is only moderately im- 
proved when varying Qbh 2 , h and Q c h 2 within the allowed 
range and the likelihood ratio also confirms that these pa- 
rameters are not very sensitive to the value of the CMB 
variance. However, this ratio reaches higher values when 
considering the variation of A, r and n s and also the oo 
becomes perfectly compatible with the model, showing that 
these three parameters are more relevant regarding this ob- 
servable. 

Taking into account the previous results, we have fur- 
ther investigated the effect of cosmological parameters by 
constructing the likelihood varying simultaneously the three 
more sensitive parameters, i.e. r, n s and A, while fixing the 
others to the values of the best-fit model. We find that there 
is a strong degeneracy between r and A, and, therefore, we 
cannot extract any further information with the observed 
CMB variance regarding these two parameters. However, 
the situation is different for the spectral index. We have 
constructed the marginalized likelihood of n s , finding that 
values of the spectral index higher than the ones given by 
the WMAP best-fit model are favoured. In particular, we 
find (n s ) = 0.98, although one should point out that, as ex- 
pected since we are using only one observable, the likelihood 
constructed in this way is significantly wider than the one 
obtained with the full power spectrum of the WMAP data. 

Therefore it is possible to find models consistent with 
the estimated CMB variance by varying the value of the 
cosmological parameters in a certain range. However, in or- 
der for a model to be valid, it should still be compatible 
with the measured power spectrum. Thus, in principle, we 
should look for the model that best fits the data, allowing to 
vary all the cosmological parameters at the same time, sub- 
ject to the constraint of compatibility with the estimated 
CMB variance. Unfortunately, in practice, this task is very 
complicated, and, as an illustration, we have just tried to 
find one model (not necessarily the best possible) that fits 
the WMAP power spectrum and is also consistent with the 
estimated CMB variance. 

In order to find this model, we have carried out the fol- 
lowing analysis. First we have selected a number (around 
20) of sets of values of {r, n B , ^4} which are contained in 
the most probable 68 per cent of the 3-dimensional likeli- 
hood previously constructed for these para meters. For each 
of these sets, we have then run CosmoMC (|Lewis fc Bridle! 
12002! ) in order to find the best-fit to the WMAP data, with 
r, n 3 and A fixed for each set, varying the rest of the cosmo- 
logical parameters. From all the considered models, we find 
that the one that best fits the WMAP power spectrum has 
the values given in Table [S] It is interesting to note that the 
value of the spectral index is higher (n s — 0.983) for the new 
model. As seen in Figure[4]the model fits reasonably well the 



Table 5. Estimated cosmological parameters for the WMAP 
best-fit and for an illustrative model more consistent with the 
estimated variance. 



Parameters 



Best-fit 



New model 



lOOf^/i 2 


9 9 o+0.07 
z - zci -0.09 


2.27 


n c h 2 


104+ 007 


0.089 


h 


n 7 n+0.03 


0.81 


T 


0.08812;™ 


0.132 


n s 


0.951™ 


0.983 


A 




0.75 




Figure 4. The power spectrum for the WMAP best-fit model 
(solid black line) and for an illustrative model more consistent 
with the estimated variance (red dashed line) are shown. For com- 
parison, the binned power spectrum measur ed from the WMAP 
three- year data (blue points) is also plotted jSgcrecl ct al.: 2007t ). 



power spectrum and the main difference can be found at low 
£'s where the new model (red dashed line) goes below the 
WMAP best-fit (black solid line) , which would tend to give 
a smaller variance. In particular, we find that the estimated 
variance has a theoretical significance of 93.1 per cent for 
this new model (to be compared with 99.3 found for the 
WMAP best-fit) and therefore the CMB variance is not so 
significantly low in this case. 

Therefore, this seems to indicate that the low value of 
the CMB dispersion found in the WMAP data could be 
made compatible with the model, while keeping a reason- 
able fit to the power spectrum, by varying the cosmological 
parameters. In addition, we would like to point out again 
that we have considered only a restricted set of models and 
if the whole set of cosmological parameters were allowed to 
vary simultaneously, it is likely that we would find mod- 
els that would fit better both the estimated variance and 
the power spectrum. Nonetheless, a more exhaustive study 
where a wider range of models were considered would be 
necessary to confirm this result. 
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5.2 Gaussianity and isotropy of the CMB 

A deviation of the CMB from Gaussianity and/or isotropy 
could also explain the inconsistency found between the esti- 
mated CMB variance and the considered model. Moreover, 
this possibility is also supported by the fact that different 
authors have claimed the presence of departures from Gaus- 
sianity and/or isotropy in the WMAP data (see 

One possible way to test the Gaussianity and isotropy 
of the CMB is to perform a consistency test of the data using 
equation Q, that relates the CMB variance and the power 
spectrum. Gaussianity and isotropy have been assumed to 
estimate both, the CMB variance and the WMAP power 
spectrum. Therefore, if these assumptions are valid, both 
results should be consistent. 

To perform this test, we have used the actual measured 
power spectrum of the WMAP dat4E i.e. the power spec- 
trum estimated for the particular realization of the CMB 
sky, instead of the best-fit model. From these CYs, we 
have calculated the power spectrum that we expect for the 
WMAP combined map (using the results of Appendix [A")) . 
Finally, we have calculated the value of the CMB variance 
from equation Q using this power spectrum of the com- 
bined map. In this way, we find a value for the dispersion 
of 8.40 x 10 -2 mK. The error associated to this quantity 
is expected to be small, since the main contribution to the 
variance is due to the low multipoles, which are almost un- 
affected by the instrumental noise . This value should be 
compared to (8.16 ± 0.04) x 10 -2 mK, the dispersion that 
we have previously estimated from the combined map. The 
difference between these two estimations indicates the pres- 
ence of an inconsistency in the data between the measured 
power spectrum and the estimated variance. 

To study further this possibility we have also calculated 
the theoretical significance of our estimated CMB variance 
from equation ([6]), using again the power spectrum of the 
WMAP combined map constructed from the actual mea- 
sured WMAP power spectrum (instead of the best-fit model 
as done in fj3]). In this case, we find a value of 96.0 per cent 
for the significance. Moreover, if we remove the quadrupole 
from the data (as explained in ^3|, the significance increases 
to 99.5 per cenfl This points out again towards an incon- 
sistency between the estimated values of the variance and 
those of the power spectrum in the WMAP data (inde- 
pendently of the best-fit model). In addition, it also shows 
that the anomaly is present at multipoles higher than the 
quadrupole. 

Therefore, a deviation from Gaussianity and/or 
isotropy of the WMAP data could be the reason for the 
inconsistency found between the estimated variance and the 

6 The power spectrum measured for the three-year WMAP data 
is available at the LAMBDA web site. 

7 In [|3] when removing the quadrupole, the significance of the 
variance decreased (unlike in the current analysis). This makes 
sense, since the quadrupole of the data is known to be signifi- 
cantly lower than that of the best-fit model, and therefore remov- 
ing the quadrupole from the analysis would make the variance 
more consistent with the best-fit model. However, in the current 
analysis, we are comparing the estimated CMB variance with the 
actual measured power spectrum of the WMAP data, where the 
quadrupole is also low and, therefore, removing the quadrupole, 
does not necessarily increase the consistency. 



considered model. Moreover, taking into account the results 
of <J4] that show that the anomalously low variance cannot 
be explained by known systematics or foregrounds (and as- 
suming that the estimated CMB power spectrum is neither 
affected by spurious signals), a cosmological origin for this 
deviation from Gaussianity and/or isotropy cannot be dis- 
carded. 



6 CONCLUSIONS 

We have studied the normalised temperature distribution of 
the CMB using the WMAP data, finding that the skewness, 
kurtosis and the Kolmogorov-Smirnov test are consistent 
with Gaussianity. However, the dispersion was significantly 
lower (at the level of 97.8 per cent) than expected from 
CMB simulations using the WMAP best-fit power spectrum. 
This result was even more significant when only the north- 
ern ecliptic sky was considered (at 99.4 per cent). In oder to 
clarify this point, a more detailed analysis of the CMB dis- 
persion of the WMAP data has been performed. In particu- 
lar, we have estimated this dispersion to be 8.16 x 10~ 2 mK, 
which is again anomalously low (98.7 per cent of the Gaus- 
sian CMB simulations have a larger value of the dispersion) . 
If only the northern ecliptic hemisphere is considered, the 
significance increases to 99.8 per cent, whereas the southern 
ecliptic hemisphere is compatible with the simulations. 

We have repeated our analysis using single radiometer 
and single year maps, finding that the significance of the 
low variance is practically frequency independent and that 
it cannot be explained by an anomalous radiometer or by 
an artifact in a single year data. In addition, the presence 
of residual foregrounds or 1/f noise cannot explain the low 
value of the CMB dispersion, since they would tend to shift 
ao towards higher values. 

We have also investigated the effect of varying the cos- 
mological parameters in the significance of the estimated 
variance. In particular, searching within a restricted number 
of models, we have fround one case that provides a reason- 
able fit to the WMAP power spectrum and for which the 
signficance of the estimated variance decreases to 93.0 per 
cent. This seems to indicate that is possible to find a model 
that is consistent with both the estimated variance and the 
measured power spectrum of the WMAP data, especially 
if a wider range of cosmological models is considered. Nev- 
ertheless, a more exhaustive study is necessary to confirm 
this result. It is also interesting to point out that the CMB 
variance tends to favour higher values of the spectral index 
than the full-power spectrum analysis of the WMAP data. 

Finally, we have considered the possibility that a devi- 
ation from Gaussianity and/or isotropy could be the reason 
for the anomalous variance. We find that there is an in- 
consistency between the estimated CMB variance and the 
actual measured power spectrum of the WMAP data, inde- 
pendently of the best-fit model. This indicates that a devi- 
ation from Gaussianity and/or isotropy could be a possible 
explanation for the found anomaly. Moreover, taking into 
account that known systematics or foregrounds cannot ex- 
plain the anomalously low CMB variance, a cosmological 
origin for this possible deviation cannot be discarded. 
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APPENDIX A: THEORETICAL POWER 
SPECTRUM OF A COMBINED MAP 

The derivation of the angular power specturm of the CMB 
signal present in the combined Q-V-W map C\ is given in 
this appendix. This combined map, T c (x), is the linear com- 
bination (see equation [TJ of the eight cleaned maps (Ti (x) 
with i running from 3 to 10) at different WMAP frequencies 
(Ql and Q2 at 41 GHz, VI and V2 at 61 GHz, and Wl, 
W2, W3 and W4 at 94 GHz). Let us recall here the explicit 
formula: 



i=8 



(Al) 



where Ti(x) is the CMB observation made by the i radiome- 
ter and the coefficients Wi(x) represent, at a given postion 
x, the relative weight given to the noise dispersion for each 
radiometer, and are normalized to unity (see equation [2} . 

Each CMB observation made by the i radiometer can 
be expressed in terms of the spherical harmonic coefficients 
of the pure CMB signal ai m , the beam (b\) and pixel (pi) 
window functions: 



Ti(x) 



l,m 



aimb\p(Yi m (x) 



(A2) 



Let us define Ti(x) — Wi(x)Ti(x) as the noise- weighted 
CMB signal observed by the radiometer i. It is straight- 
forward to show that the angular power spectrum of the 
combined CMB signal is given by: 



Cf 



i 



\ ^ 



21 + 1 ^ 



(A3) 



where C\ 3 is the cross-angular power spectrum of the CMB 
signal between the i and the j radiometers and a\ m axe the 
spherical harmonic coeficients of the noise-weighted CMB 
signal observed by the radiometer i: 

dQ fi(x)Yl m (x) = J dQ Wi(x)Ti(S)Yl m {x) = 

= E <V ml K' Pi' J dn w i ( f ) Y im (2) Y t > m , (x) = 
e'm' 

= y^o,/ ibitP./JCt „/ i. (A4) 



K l , i is a matrix introducing coupling among multipoles 
due to the noise- weighting. 



Including equation (|A4[I in equation (|A3|I . the the an- 
gular power spectrum of the combined CMB signal reads: 



a\ m a J l71 



C * = EaZTlE' 

i,j m 

= E27TlE E a t',m' b l'Pi' K trra' m < 
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(A5) 



where C,i is the angular power spectrum of the pure CMB 
signal and we have applied TV/ / // a„i ia*,< ,, — 

E„i 1 C „i 5„i „" 5 1 a . 
i ,m £ t i mm 



The coupling K) 



can be expressed in terms of the 



Wigner 3j-symbols and the spherical harmonic coefficients 
(m( ra ) of the noise-weights, as follows: 



K 



; iV = / dQ Wi{x)Y; m (x)Y t , m ,(x) = 

f dQ F/ ro (£)F /m ,(f)l> m ,,(s) 



/ j £ m 

l" ,m" 




One of the ortho gonality properties of the Wigner-3j 
symbols reads (see e.g. iHivon et al.ll2002l ): 
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m l j rn 2 
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. >l ., + i 5t 3 e i 5 m3mi 6(h,t2,£ 3 ) (A7) 

where 6(li, h, h) = 1 when the triangular relation \£i —£z\ ^ 
£3 ^ £1 + £2 is satisfied. 

Replacing equation (|A6fl in equation (| A5|) and taking 
into account equation (|A7I) . it is straightforward to derive 
the final expression for the angular power spectrum of the 
combined CMB signal: 



ci = EE^' 6 ^ 

i,j l' 



21+1 
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47T 

£ £ e" 




(A8) 



where w % j, is the cross-angular power spectrum of the noise- 
weights between the radiometers i and j. 

In order to test the performance of this approach, the 
theoretical angular power spectrum of the combined Q-V- 
W WMAP three-year map (at n S ide=256 resolution) given 
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Figure Al. Relative error (solid black line) of the theoretical 
angular power spectrum of the WMAP combined map, given by 
equation ]A8t . with respect to the one averaged from 1000 sim- 
ulations. The relative error is obtained as the difference of the 
theoretical Cg's minus the average ones, and this difference is di- 
vided by the average C/s. The dashed red lines give the lcr errors 
(obtained from the simulations) divided by the averaged angular 
power spectrum from simulations. 



by the above equation has been compared with the angular 
power spectrum obtained from the average of 1000 Gaussian 
simulations of the combined map. In Figure |A"T1 the relative 
error between both spectra is given (solid black line) up to 
£ = 2.5n si d e . The red dashed lines give the relative lcr er- 
ror (obtained from the simulations). The agreement between 
both angular power spectra in the considered range is very 
good. 



